www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第三讲 蒙特卡洛方法基本思想/liti2.m

    function liti2(n,p,mm)
pro=zeros(1,mm);
randnum = binornd(n,p,2,mm);a=0;
 for i=1:mm
 a=a+randnum(1,i)*randnum(2,i);   pro(i)=a/i;
end    
pro=pro,num=1:mm;plot(num,pro)